Compression of nanowires using a fiat indenter: Elasticity measurement in nanoscale 
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A new experimental approach for the characterization of the lateral elastic modulus of individual 
nanowires is demonstrated by implementing a micro/nano scale diametrical compression test geom- 
etry, using a flat punch indenter inside of a scanning electron microscope (SEM). A 250 nm diameter 
single crystal silicon nanowire is tested. Since silicon is highly anisotropic, the compression axis of 
the wire was determined by electron backscatter diffraction (EBSD). A two dimensional analytical 
closed-form solution based on a Hertz model is presented. The results of the analytical model are 
compared with those of finite-element simulations and to the experimental diametral compression 
results and show good agreement, (revised manuscript) 



Nanowires have attracted great interest in the past 
decade due to their interesting mechanical, electrical, and 
optical properties [TJ [2] . This has enabled nanowires to 
gain advantages in a wide range of applications such 
as nanoelectromechanical systems [3], gas sensors [I], 
resonators [5], light sensors 6 , photovoltaic cells [TJ, 
and atomic force microscopy and scanning tunneling mi- 
croscopy tips [5|. In this paper we focus on the charac- 
terization of the mechanical properties of nanowires, in 
particular by in-situ scanning electron microscopy (SEM) 
diametral compression using a flat punch. 

The main experimental techniques employed to de- 
termine the mechanical properties of nanowires are me- 
chanical resonant methods [§1 [TU] , AFM-based bending 
methods [TJ rTJJ [T2] , tensile testing [T3J , and nanoindenta- 
tion methods [T1HT5] . As pointed out by several research 
groups, nanowire indentation has the great advantage of 
little specimen preparation, in fact the nanowire does not 
need to be held in place as with all the other techniques. 
Not only does the holding require intensive preparation 
but it can also influence the measured mechanical re- 
sponse since the mechanical and adhesive properties of 
the holding element are generally neglected. Nanoinden- 
tation has a couple of drawbacks such as the alignment 
of the nanoindenter tip with respect to the nanowire [17] 
and the non-applicability of the Oliver and Pharr method 
[T51 HI?] since the nanowire cannot be considered a half- 
space [20]. 

To overcome these issues, it is introduced here the in- 
situ scanning electron microscope (SEM) diametral com- 
pression test performed at the mirco/nano scale. In this 
approach, the misalignment defined as the non-perfect 
perpendicularity between the nanowire axis and the nor- 
mal to the flat punch is attenuated since the experiments 
are conducted in a SEM while the misalignment due to 
the non-intersection of the nanowire axis and the flat 
punch axis is not as critical as in the case of wire inden- 
tation p~7] since the indenter used is a flat punch with 
a diameter that is much larger than the diameter of the 
nanowire. Concerning the diametrical compression data 



analysis, Hertz theory is well suited. This allows one to 
solve, at least in the 2D simplified case described in the 
next section, in closed-form this problem allowing one to 
determine the lateral elastic modulus of the wire from the 
load-displacement data obtained experimentally. In or- 
der to validate this technique, the experimental test was 
compared to an analytically closed-form and a numerical 
solution of the diametrical compression test problem. 

To the authors knowledge this is the first time that 
a wire is compressed by a flat punch in the micro/nano 
scale. This test geometry is analogus to the macro-scale 
Brazilian test [3TJ [23] that was introduced back in 1950s 
to study ceramic fracture behavior. 

The remainder of this paper is organized as follows. 
In Sec. II we first describe the formulation of our analyt- 
ical approach. The results of the analytical model, ex- 
periments and numerical simulations are compared and 
discussed in Sec. Ill, followed by conclusions in Sec. IV. 

In this section, we outline a closed-form Hertzian 
solution to determine the elastic modulus of the 
nanowire tested with a flat indenter, analyzing the load- 
displacement (P—S) relationship. We study the cross sec- 
tion of the system in two physical dimensions. The total 
displacement 5 measured from the compression experi- 
ment consists of three terms: S — Si + 5 W + S s where the 
subscripts i, w and s indicate the indenter, the nanowire 
and the substrate, respectively. 

To calculate the deformation, we seek the two compo- 
nents of stress (a x ,a y ) distributed at an arbitrary point 
A(x, y) at x = in the cylinder (of radius R) under a 
pressure P (Figure[T]). By assuming a Hertzian traction 
distribution (= 2Py/l — x 2 /a 2 /air and R >> a) at the 
interface [531, the stress can be written as 
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where a is the transverse semi-contact length (Figure]!]), 
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which can be calculated according to a Hertzian distri- 
bution of pressure at the interface as 



In summary, the load-displacement relationship can be 
expressed as 



y/APR/TTE i 



(2) 



where E mt stands for the interface elastic modulus, 
which is defined as 



= 1/E* + 1/E* 
l/EZ* = l/E* s + 1/E, 



(3) 



where the subscripts iw and sw indicate the indenter- 
wire and the wire-substrate interfaces, respectively. E* = 
E/(l — v 2 ) where E is the elastic modulus and v is the 
Poisson's ratio. 
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FIG. 1: 2D diagram of an indentation test on a nanowire 
(cylindrical cross section). 

The vertical deformation e y at A(x = 0, y) under biax- 
ial stress is 
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(4) 



The compression of the wire S w is given by integrat- 
ing e y in each half cylinder e y dy. The result of the 
integration (using Eqs.(l-5)) gives 
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Additionally, the penetration of the nanowire into the 
indenter and the substrate refers to a typical problem 
of calculating the compression of a half-space under a 
Hertzian pressure, assuming that the dimension of the 
substrate and the indenter is much larger than that of the 
contact interface. A classical solution [23 to this problem 
involves integration along the depth of the substrate (h s ) 
or indenter (hi). It yields: 
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where h is thickness of the indenter or the substrate. 
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It is worth noting that this relationship between the 
load P and the displacement 5 is not totally linear since 
a depends on P. 

The compression of nanowires is simulated using a 2D 
finite-element (FE) model realized with a commercial 
computation code [53]. The geometry of the simulated 
system is shown in Figure]!] right. The bottom of the 
substrate is fixed in the simulations and an imposed ver- 
tical force (downward) is applied to the top of the punch 
indenter. The object sizes and the material properties of 
the simulation system are set with respect to our exper- 
imental setup in nanoscale. For the mesh, a plane-strain 
Lagrangian model is used (approxiamatly 40000 triangle 
meshes). These interface model is set to be the general 
contact model without friction. 

The 250 nm diameter Si nanowire used in our experi- 
ments was prepared by the VLS method using chemical 
vapor deposition (CVD) with silane. The substrate was 
heated to about 580° C at a pressure of 2 x 10~ 7 mbar 
for 10 minutes. The temperature was then reduced to 
510°C and a mixture of Ar (10 seem) and S1H4 (5 seem) 
was introduced for the nanowire growth for 15 minutes 
at a pressure of 2 mbar. 

The nanowires were compressed with a Hysitron 
Picolndenter® that was built for SEM in-situ experi- 
ments [25 . The experiments were conducted in a Zeiss 
DSM 962 SEM operated at an accelerating voltage of 
10.0 kV. All compressions were conducted such that the 
substrate tilt was between 8° to 10° with respect to the 
electron beam. The tip was a boron-doped diamond that 
was focus ion beam (FIB) milled to create a flat-punch ge- 
ometry. The nominal indenter tip displacement rate was 
2 nm/s using a closed feedback- loop which adequately 
represents displacement control during elastic segment 
of the compressions. An SEM image of the experimental 
configuration can be seen in Fig. 2(a). 

The analytical model (Eq.[7]) for the deformation of the 
nanowire requires accurate values of the Young's Modu- 
lus E and the Poisson coefficient v for all involved materi- 
als. However, the elastic behavior of single crystal silicon 
is anisotropic and strongly depends on the orientation of 
the crystal [26\ . In order to determine the effective values 
of E and v in our experiments, the crystallographic ori- 
entation of the compression axis of the specific nanowire 
has been determined by means of electron backscatter 
diffraction (EBSD), which can be used to identify the lo- 
cal crystal orientation with a spatial resolution as good as 
10 nm and an angular orientation resolution better than 
0.5° [37]. The EBSD analysis was conducted on a Hitachi 
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FIG. 2: (Color online) (a) SEM image showing a nanowire 
lying on the Si wafer substrate and the diamond flat-punch 
indenter tip. (b) Streographic projection of the <100> di- 
rections illustrating the crystal orientation of the Si nanowire 
(blue dots) and the substrate (red dots), (c) EBSD mapping 
of the Si nanowire and the substrate. The color code gives 
the crystallographic orientation of the surface normal. 

S-4800 cold field emission high resolution SEM equipped 
with a TSL/EDAX EBSD system at an acceleration volt- 
age of 20.0 kV. The crystal orientation mappings of the 
Si nanowire and the underlying Si wafer can be seen in 
Figure[2} 

Two consecutive compressions were conducted on the 
wire. The first compression was elastic such that there 
was minimal difference between the loading and unload- 
ing slopes of the load-displacement graph seen in Fig- 
ure[3] A second compression was then conducted at 
which point plasticity occurred. The plasticity enabled 
the length of the contact to be estimated as 1.34 /zm by 
subsequent SEM plan-view imaging of the deformed re- 
gion of the nanowire. This longitudinal contact length 
is very important since it must be used to normalize the 
load, P, in Eq.[7] It should be noted that the contact 
length within the first 10 nm of displacement in Figure[3] 
is not constant, but increases to the final recorded length 
of 1.34 fim. This is due to small but unavoidable mis- 
alignments between the flat punch, nanowire and sub- 
strate. Excluding these first few nanometers, the remain- 
der of the elastic portion of the first compression is then 
used in order to fit the effective elastic modulus of the 
nanowire using Eq.[7] 

The EBSD mapping confirms that our Si nanowire is 
a single crystal. The crystal orientation of the nanowire 
and the underlying substrate are visualized in the pole 
figure of Figure[2j The <100> crystal directions of the 
nanowire (red) and the Si wafer substrate (blue) are dis- 
played in a stereographic projection in Figure[2^b) . The 
configuration of the nanowire is drawn schematically into 
the pole figure. The crystal orientation of the nanowire 
is such that the compression direction y expressed in the 
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FIG. 3: (Color online) Comparison of the load-displacement 
curves from a nanoindentation experiment (load: Elastoplas- 
tic and unload: Elastic), the analytical model (Eq.[7| and a 
FE simulation. 

crystal coordinate system is (0.17, 0.92, 0.34). The lat- 
eral direction x is (0.87, -0.39, 0.29). According to the 
measured crystal orientations, the Young's modulus and 
the Poisson's ratio of the nanowire are estimated to be 
E w — 147 GPa and tv,=0.233 based on the data of Wort- 
man and Evans |26) . In the substrate, the principal di- 
rection of the stress varies and strongly depends on the 
position within the substrate. Therefore, based on in- 
dentation results, an average value of 153 GPa has been 
assumed for the modulus of the substrate [28] . 

Figure[4] shows a comparison between the compres- 
sion experiment, the theoretical prediction and the finite- 
element (FE) simulation. It can be seen that the experi- 
mental elastic unloading curve roughly follows the track 
predicted by the simulation using the effective Young's 
modulus and Poisson's ratio estimated with EBSD. At 
the small deformation range, the load-displacement ra- 
tio predicted by the Hertzian model is close to that of 
the experiment and simulation. The difference becomes 
more important when the deformation increases. 

The stability of this simplified Hertzian model has been 
checked by using FE simulations [24 . Our experimental 
configuration (flat punch + wire + substrate) was rebuilt 
in a plane-strain Lagrangian model applied to 4 different 
sets of material parameters (see Table 1). A vertical load 
was applied progressively to the top of the indenter (see 
FigureJT] right) and the deformation was then measured 
to compare with that predicted by Eq.[7] The comparison 
of the P — 8 curves is shown in Figure^] We can see that 
the stiffness increases with increasing elastic modulus of 
the materials. The stiffness obtained by simulations is 
sightly lower than that predicted analytically when the 
deformation is small. However, the difference becomes 
more significant in the case of large deformation. This 
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FIG. 4: (Color online) (a) Left: Geometry meshes used in our 
FE indentation simulations. Color scale represents the distri- 
bution of strain energy (N-m/m). A real-time video-recording 
of simulation is available online 29 . (b) Comparison between 
the load-displacement curve from Eq.[7] (lines) and that from 
simulations (symbols) with 4 different sets of parameters as 
that shown in the inset: Young's modulus E (GPa) and Pois- 
son's ratio v {hi = h s = lOOi? = 10 /im). 



may be due to the fact that the R >> a Hertzian assump- 
tion is no longer valid in the large deformation case. 

In conclusion, in-situ SEM diametrical compression at 
the nano-scale has been introduced to measure the elastic 
modulus of individual nanowires with cylindrical cross- 
sections. The experimental diametral compression load- 
displacement response of a 250 nm diameter single crys- 
tal Si nanowire was compared to an analytical solution. 
The analytical solution of the three element contact me- 
chanics problem was derived from Hertzian theory where 
the geometry of the experiment is simplified as the com- 
pression of an infinite long cylinder (nanowire) between 
two infinite half-spaces, the diamond flat punch and the 
silicon wafer respectively. Well characterized materials 
(diamond, silicon wafer and nanowire) were chosen for 
the three contact elements in order to have known ma- 
terial property values to feed into the analytical model. 
The comparison between the experimental measurement 
and the analytical solution showed good agreement and 
thus the feasibility of such experiment to determine the 



lateral elastic modulus of nanowires. The stability and 
reliability of the analytical solution was proven by com- 
paring the load-displacement response obtained by the 
analytical model with the response obtained by FE nu- 
merical simulation applied to four different set of material 
parameters. 
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